Markarian 6: shocking the environment of an intermediate Seyfert 
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ABSTRACT 

Markarian 6 is a nearby (D~78 Mpc) Seyfert 1.5, early-type galaxy, with a double set of 
radio bubbles. The outer set spans ~7.5 kpc and is expanding into the halo regions of the host 
galaxy. We present an analysis of our new Chandra observation, together with archival XMM- 
Newton data, to look for evidence of emission from shocked gas around the external radio 
bubbles, both from spatially resolved regions in Chandra and from spectral analysis of the 
XMM data. We also look for evidence of a variable absorbing column along our line of sight 
to Mrk 6, to explain the evident differences seen in the AGN spectra from the various, non- 
contemporaneous, observations. We find that the variable absorption hypothesis explains the 
differences between the Chandra and XMM spectra, with the Chandra spectrum being heavily 
absorbed. The intrinsic Nh varies from - 8x10^' atoms cm~^ to ~ 3x10^^ atoms cm"^ on short 
timescales (2-6 years). The past evolution of the source suggests this is probably caused by a 
clump of gas close to the central AGN, passing in front of us at the moment of the observation. 
Shells of thermal X-ray emission are detected around the radio bubbles, with a temperature of 
~0.9 keV. We estimate a temperature of ~0.2 keV for the external medium using luminosity 
constraints from our Chandra image. We analyse these results using the Rankine-Hugoniot 
shock jump conditions, and obtain a Mach number of ~3.9, compatible with a scenario in which 
the gas in the shells is inducing a strong shock in the surrounding ISM. This could be the third 
clear detection of strong shocks produced by a radio-powerful Seyfert galaxy. These results 
are compatible with previous findings on Centaurus A and NGC 3801, supporting a picture in 
which these AGN-driven outflows play an important role in the environment and evolution of 
the host galaxy. 

Subject headings: Active Galaxies - Radio Galaxies - Mrk 6 - X-rays - Chandra - XMM- 
Newton - 
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1. Introduction 



Recent Chandra o bserva tions of the environ ments of sev eral powerful radio galaxies (e.g. Hydra A, 



McNamara et al 



McNamara & Nulsen 



2OOOI: M8 7. lYoung etal.ll2002l : Hercules A. iNulsen et all 120051 : see also the review by 



20071) have led to significant progress in understandi ng the AGN-driven gas outflows 



in the se systems and the role they play in galaxy formation and evolution dCroton et al.ll2006l : 
20o3). 



Bower et al. 



We now know that, although the most powerful radio outflows, spanning hundreds of kpc, are associ- 
ated with massive elliptical systems, smaller structures also connected to an active n ucleus can be found in 
a variety of systems and en vironments , including spiral and di sk galaxies (see e.g. IGallimore et all 12006 , 
Kharb et al.ll2006l.lHota & Saikia.2006l . lSaikia & Jamrozyll2009l) . The mechanism by which these structures 
are produced is likely to be related to the one we see in the most powerful sources, but on a smaller scale. 
Most of the observed AGN-driven bubbles have been found to be overpressured with respect to their sur- 
roundings, and may be thus inducing shocks into their surrounding medium. The wide range of morpholo- 
gies of the galaxies wh ere radio bubbles have b een found, and the fact that this AGN-driven phenomenon 
is most likely episodic (|Saikia & Jamrozyll2009r ) make understanding the energetics involved in this process 
fundamental to estimate its impact on AGN feedback and galaxy evolution, and extrapolate how common 
this mechanism can be in low-power systems. 

In our search for evidence of galaxy feedback associated with kpc-scale radio bubbles, we have carried 
out obse rvations of a variety of systems, the most no table perhaps being the nearby Fanar off-Riley type 
I (FR I, Ipanaroff & Rileylll974) galaxy Centaurus A (|Kraft et all boosl fcroston et allboogl). finding evi- 



dence for shocks also in smaller, more distant systems su ch as NGC 3801 (ICrqston et al. 



20071) and a rather 



more complex scenario in the spiral galaxy NGC 6764 (|Croston et allbOOSbl ). We recently carried out a 
Chandra observation of M arkarian 6 (Mrk 6, IC 450, z = 0.018676), an early-type SO Seyfert 1.5 galaxy 
(lOsterbrock & Koskilll976h whose characteristics have been studied in the radio iLi 4GHz = 1-7 x 10^^ W 
Hz~^ sr~'), infrared, optical and X-ray wavelengths over the last 30 years. Recent radio stu dies have un 



veiled a complex structure surrounding the AGN, with a doubl e set of bubbles and a radio jet (IKukula et al. 



19961) . suggesting a jet precession scenario (IKharb et al.ll2006r) . In this paper we describe the results of our 
analysis of the Chandra observation. 

We have also analysed three previous datasets from XMM-Newton to study the evolutio n of the AGN 



and i t s immediate enviro nment over time. These data were analysed in detail previously (|Schurch et al 



20061 : llmmler et all 120031) . using difi'erent models to address the complex scenario surrounding the AGN. 
We have approached this analysis by searching for consistency between the Cha ndra and XMM-Newton 
datasets in the context of variable absorption, following both the methods used by iHardcastle et al.l (120091) 
when modelling a sample of 3CRR radio sources, and a double partial covering model that has previously 
been successful in describing the properties of the nuclear spectrum of Mrk 6. 



Throughout this paper we use a concordance cosmology, with Hq-10 km s ^ Mpc \ n„ 
nA=0.7. The redshift of Mrk 6 corresponds to a luminosity distance of 78.4 Mpc. 



-0.3 and 



-3- 



2. Observations and data reduction 
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Fig. 1. — Left: Chandra A CIS-S 0.4-3 keV im age of Mrk 6, with cr=3 pixels Gaussian smoothing and 
overlaid VLA radio contours (iKharb et al. showing the match between the radio and X-ray emission. 

In this frame are also displayed the source (blue, continuous) and background (magenta, dashed) extraction 
regions for the bubble spectra. The label A shows a possible background quasar, while B indicates the 
emission coincident with the extended [OIII] emission-line regions, see Figure |2l The energy range has 
been selected to highlight the shell-like structure surrounding the bubbles and the overlap between the X-ray 
emission and the radio structure. Right: (non-smoothed) XMM-Newton's MOS 1 image (obsid 0305600501) 
of Mrk 6, with the source (blue, continuous) and background (magenta, dashed) extraction regions for all 
the AGN spectra {Chandra and XMM). 



We observed Mrk 6 with Chandra's ACIS-S CCD on July 12 2009, with a total exposure of 75 ksec. 
The observation was carried out in faint mode, to minimise issues with the background and at the same time 
avoid telemetry saturation. There were no intervals of high background during this observation, hence we 
analyse the full exposure. We reprocessed the data from level 1 events in the standard manner, using CIAO 
4.2 and CALDB 4.2.0. To improve the quality of the data for the purposes of a study of the extended radio 
bubbles, we removed the pixel randomization and ran the destreaking routines on the data. We produced 
two filtered images (0.3-7 keV and 0.5-5 keV) to assess the extent and structure of the radio bubbles, and 
chose our regions accordingly, but extracted all the spectra from the original events file. We performed our 
spectral fits with XSPEC, constraining the energy range to be coincident with that covered by the instrument 
calibration (0.3-7 keV). We estimated the pileup fraction for these data from the Chandra documentation 
and calculated a numerical value using the PIMMS tools. The pileup fraction is about 12% throughout the 
observation (there are no substantial differences in the count rate during the observation). 

The XMM-Newton observations were taken in March 2001 (obsids 0061540101 and 0061540201), 
April 2003 (obsid 0144230101) and October 2005 (obsid 0305600501). We reduced the data using the 
standard routines from SAS version 9.0 and the latest calibration files. Table[T]gives details on the exposure 
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times. The 2001 and 2003 observations were taken with the medium filter, while 0305600501 was taken with 
the thin filter. We decided not to use the observation 0061540201 because it was taken in calibration closed 
mode. We discarded the first 10 ksec from the observation 0144230101 due to high background. Despite the 
background contribution being uneven during this observation, after examining the relative count rates of 
the source and the background, we decided not to discard any other time intervals, since the background in 
these periods was always below 10% of the intensity of the source across the whole energy range, which is 
accurate enough for our purposes. We extracted spectra for the PN, MOS 1 and M0S2 cameras, except in the 
case of the observation 0144230101, where a PN spectrum could not be obtained. We limited our spectral 
fits to the 0.3-8.0 keV energy range, to analyse an energy range comparable to that covered by Chandra. 

For the study of the properties of the AGN we used the same extraction region for the Chandra and 
XMM data, a ~60 arcsec radius circle centred in the source, which contains most of the MOS and PN PSF. 
Although we know from the Chandra image that there are a few point sources other than the AGN within 
this extraction region (see Figure [T|), they are not resolved by the XMM instruments, and their intensities are 
so low compared to the AGN that we can consider their effects negligible. We also use the same background 
region for all the data, a -50 arcsec radius circle. North of the source region to avoid contamination from 
the host galaxy, as seen on Figure [T] 



To simulate the AGN on the Chandra detector we used the ChaRT 1.0 web interface and MARX version 



4.5. 




Fig. 2. — [OIII] image of Mrk 6 from iKukula et al.l (119961) . with overlaid radio contours, showing the 
structure of the extended emission-Une region (ENLR). Scale: 1"=366 pc. 



For the study of the extended bubble area, we used as a reference the radio data from iKharb et al. 



(l2006h . defined contours around the largest-scale bubbles and drew source regions outside these contours, 
keeping well away from the central AGN to avoid contamination from its PSF, as well as from the inner 
radio bubbles and jet. Figure [T] shows details on the extraction regions and the structure of the extended 
emission from the bubbles. In this case, we took care to exclude the only likely background source (South 



-5- 



of the source, see discussion in Section [3T2b in order to avoid contamination in our spectra. 



All tlie spectra were grouped to 20 counts per bin after background subtraction, before spectral fitting, 
to be consistent with st atistics. We used a fixed Galactic absorption of Nh = 6.39 x lO^^'atoms cm~^ 
(jpickey & Lockman|[l99ol) for all our spectral fits, and a redshift z ^ 0.018676 (SIMBAD). 



Table 1: Summary of Mrk 6 X-ray observations. 



Telescope 


Instruments 


Date 


Obsid 


ksec 


XMM-Newton 


PN, MOSl,2 


2001/04/26 


0061540101 


31 


XMM-Newton 


M0S1,2 


2003/04/26 


0144230101 


41 


XMM-Newton 


PN, MOSl,2 


2005/10/27 


0305600501 


20 


Chandra 


ACIS 


2009/06/12 


10324 


75 



All the parameters estimated with XSPEC are quoted with 90% confidence uncertainties. The errors 
on the upper limit of counts for the extended regions are Icr. 



3. Results 



3.1. The Seyfert core 

3.1.1. Instrumental corrections 

When we first approached the analysis of our data, we extracted a spectrum using the regions illustrated 
in the right panel of Figure [T] and fitted a model consisting of a double power law with Galactic and heavy 
intrinsic absorption. We were surprised to find that the Chandra and XMM-Newton spectra looked radically 
different. After checking the observation details of the Chandra data, using the on-line tool PIMMS , we 
found that the AGN spectrum had a noticeable pileup fraction (~12%), compatible with the frame time of 
the observation (3.1 s) and the background-subtracted count rate of the AGN (~0.11 counts/s). This resulted 
in a "ghost" peak around 2 keV in our spectra, caused by piled up 1 keV photons. We could not obtain 
satisfactory results from the XSPEC pileup model, perhaps due to the complexity of our underlying model, 
so we decided to attem pt a different approach , previously applied with success by iGetman et al.l (l2005h . 



Evans et al 



(120051) . and lHardcastle et all (120061) . 



In most sources only the central pixels of the PSF are affected by pileup. These innermost pixels also 
contain a large fraction of the total photon counts; so eliminating them effectively removes any issues with 
pileup, but it also results in a lower signal to noise ratio. If pileup is significant enough to cause a bias in the 
spectral fits, however, removing its effects often results in improved fits. In our case, excluding the central 
four pixels of the PSF effectively eradicates any pileup issues in our spectra; we extracted a spectrum from 
such a region, centered on the core, and generated an RMF and ARE in the standard way. However, since the 
PSE is energy-dependent, excluding these pixels and extracting a spectrum from the resulting annular region 
alone does not solve the problem; it is necessary to correct the ARE to account for this energy dependence. 
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Fig. 3. — Left: real Chandra ACIS-S (75 ks) 0.3-7 keV image of Mrk 6, with cr=2 pixels Gaussian smooth- 
ing. Right: ChaRT/MARX simulated data (150 ks) 0.3-7 keV, with o"=2 pixels Gaussian smoothing. The 
scales and colour bars in both pictures are the same. The simulated image contains only AGN photons, 
with no background, and helps us correct the pileup effects in our AGN spectra (we did not simulate pileup) 
and establish a statistical contribution of photons from the AGN to the bubbles and external regions of the 
galaxy. Since the exposure time of the simulation is twice as long as the one of the original data, we have 
corrected for any exposure-dependent effects when using it for our analysis. The equivalent of the faint 
streak in the simulated image is below the background level in the real data. Scale: 1"=366 pc. 



To correct the ARF we followed the procedure described by lHardcastle et al.l (|2006r) . We generated an 
energy (keV) versus flux (photons/cm^/s) table from our model fit and fed it to ChaRT, the Chandra Ray 
Tracer dCarter et al.ll2003l) . a code that simulates a PSF from a specified model. We selected an exposure 
time of 150 ks, twice as long as our real exposure, to get the best possible photon statistics in the simulated 
data while trying not to exceed the limit on the ray density. We did not simulate the pileup with ChaRT. The 
output from ChaRT is a set of rays that cannot be directly used for analysis: they need to be projected on the 
detector plane and converted into an events file using MARxQ. 

We then generated a new events file from our original data and an annular extraction region, identical 
to the one we used to generate our spectra, but excluding the central four pixels. We used a code to fit a 
5-th degree polynomial to the ratio of this events file and the whole simulated events file as a function of 
energy. This code reads in the ARF generated by CIAO and scales the effective area at each energy, using 
the polynomial fit, to effectively correct for the missing effective area due to the exclusion of the central 
pixels. The code then writes a new ARF which can be used to correct for the effects of excluding the central 
pixels. 



'See |http://space.mit.edu/CXCA^RX/docs.html] 
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We were then able to carry out fits to our extracted spectrum. As expected, the loss of counts from 
the central pixels slightly decreased the signal to noise ratio, but eliminating the bias caused by pileup 
resulted in a noticeably improved fit (see Section [3.1.2l for details) and we obtained a better statistical result. 
The simulated image is displayed on the right panel of Figure [3l next to the original data. Although the 
simulation predicts some streaking in the image, in our real exposure it is so weak it is not detectable over 
the background noise, so that we cannot use it to constrain the AGN spectrum. 



3.1.2. Models 



Early AS CA data on Mrk 6 have already made clear the difficulties underlying the study of intermediate 
type Seyferts (IFeldmeier et al.lll999l) . particularly when disentangling the absorption component from the 
intrinsic continuum shape, and several attempts have been made since then to successfully model both 
the underlying, complex physical scenario, and the properties of the X-ray emission we see. The partial 
co vering model used fo r the ASCA data was also employed in the analysis of an e arly XMM-Newton dataset 
by llmmler et all (l2003b . while analysis of a later observation (iSchurch et al.ll2006h favoured the inclusion of 
a reflection component. 

After correcting for pileup in our Chandra spectrum, it still looked very different to the XMM-Newton 
spectra. Moreover, we encountered some problems when fitting the Chandra spectrum alone, using a double 
power law model with local and intrinsic absorption (see below for a detailed description of the model). 
XSPEC is unable to disentangle the contribution to the model from the power law from that of the absorbing 
column, which results in extremely low or even negative values of the photon index for the second power 
law. Seeing that the Chandra spectrum has fewer relative counts and a different shape in the 2-5 keV range, 
we decided to test the hypothesis of a variable absorption column along our line of sight to Mrk 6, which may 
be caused by the movement of clumpy gas from the regions within a few pc of the black hole, happening on 
timescales of months t o years. This hypot hesis has been successfully used over a variety of S eyfert 1 and 2 
X-ray spectra (see e.g. iRisaliti et al. Il2002h and has already been suggested before for Mrk 6 by llmmler et all 



(|2003[) after they observed a substantial change in the absorption column between their XMM and BeppoSAX 
observations. Ho wever, the absorp tion variations between the XMM spectra from 2001, 2003 and 2005 are 
much smaller (Schurch et al.ll2006l ). Disentangling intrinsic AGN variability from the variable obscuration 
is a difficult task, both in the X-rays and at other wavelengths. We know from the optical wavelengths that 
changes in intensity of the Ha and HyS optical emission lines are a good probe of the variability of the central 
AGN itself, but there is also evi dence that suggests tha t the gas where these lines are produced can undergo 
substantial variations over time (IRosenblatt et al.ll 19921) . 



We d ecided to approach th e modeling of the AGN both from the perspective used in iHardcastle et al. 



(120091) and lCroston et al.l(|2007l) . which has proved successful in describing the properties of many powerful 
radio sources at low to intermediate redshifts, and under the parti al covering mod els that have been success- 
ful at describing not only the properties of the spectrum of Mrk 6 (|Feldmeier et al..J999, : iImmler et al ] l2003l : 
Schurch et al.ll2006l) . but also several other classical Seyferts. 
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Fig. 4.— Energy-unfolded XMM-Newton PN (black: 0061540101; red: 0305600501) and Chandra ACIS-S 
(green) spectra of Mrk 6 to illustrate the changing spectral shape caused by the variation in N//. The model 
is a power law with local absorption, with the photon index and normalization fixed to the values obtained 
in the fit of the second model (see Table |2]). 



The first model provides an accurate description of the nuclear spectrum of Cen A (|Evans et alJl2004h 
and is often used to fit spectra of narrow-lined radio galaxies. We fitted a single power law with a fixed 
Galactic absorption to the data {wabs XSPEC model). We added to this model a second absorption column 
{zwabs XSPEC model) at the source's redshift, which was allowed to vary, and a second power law. We 
added a redshifted Gaussian to account for the Fe Ka emission line. See Figured for an overview of the 
Chandra spectrum. 



Table 2: Best fitting parameters for the Chandra and XMM-Newton AGN fit. Model 1: wabs(apec+po+zwabs(po+zgauss)); Model 2: 
wabs(apec+zpcfabsi *zpcfabs2(po+zgauss)). Both models include the contribution from the thermal emission from the radio bubbles (see 
Section [X2] for details). 
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Fig. 5. — Chandra AGN spectrum with the best fitting parameters for Model 1 (see Table|2]for details). The 
model in red does not take into account the thermal emission from the radio bubbles, the model in blue does, 
and fits the data better around 1 keV. The addition of the thermal component increases the model predicted 
flux by 12% in the 0.7-1.1 keV range. 

Although the Chandra data have the spatial resolution required to differentiate the AGN contribution 
from the extended thermal emission of the radio bubbles, we used the same extraction regions on all our 
spectra (see Figure [T]), so as to be able to make a direct comparison between them. We decided to take into 
account the contribution of the thermal emission from the bubbles in our spectra by adding the apec XSPEC 
component to our model. We will discuss the details and implications of this extended emission in Section 
13.21 an overview of its effect on our AGN spectrum can be seen in Figured 

We individually fitted all the spectra, for each instrument of the XMM and Chandra, and then tied 
together the parameters for the XMM spectra belonging to the same observation, to assess the number of 
variable parameters: we assumed that the normalizations do not vary within the same observation between 
the MOS and PN instruments, and tied together the values for the three instruments, with no significant 
impact on the final statistical result. We finally attempted a full joint fit with four groups of data, one for 
Chandra and one for each XMM observation. When setting the variable parameters for the joint fit, we also 
decided that the photon indices for the second power law should be tied across the datasets, to break possible 
degeneracies with the local absorption columns, and froze the first photon index to a reasonable value (1.5) 
to minimise the overlap degeneracy between both power laws. This allowed us to make the joint {Chandra 
+ XMM) fitting considerably faster. See Figures [5] and [6] for an overview of the Chandra and XMM spectra 
and Table[2]for details on the fit parameters for this model (Model 1). 

The initial fit to the data resulted in;^'^=3331 for 3149 degrees of freedom (reduced ;^'^= 1.058). After 
correcting for the pileup effects, with the methods described in Section [3.1.11 we obtained ^^=3125 for 
2980 degrees of freedom (reduced 1.049). After adding the contribution from the radio bubbles, with 
the parameters frozen to their best fit values, we obtained =3 116 for 2980 degrees of freedom (reduced 

1.046). Since the corrections mainly affect the Chandra spectrum their contribution to the global result 
is quite small, but not negligible. 
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Fig. 6. — XMM AGN spectra with the best fitting parameters for Model 1 (two power laws + intrinsic 
absorption + thermal emission from the bubbles, see Table |2] for details). From top to bottom, obsids 
0061540101, 0144230101, 0305600501, with the PN spectrum in black, MOSl in red and M0S2 in green. 

As shown in Tabled we obtained a combined value for the index of the second power law of 1.16;J;q q3, 
which is quite flat. This is partly caused by the difficulty of disentangling the index of the power law from 
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the heavy absorption in the Chandra spectrum, as mentioned above. We find and fit the Fe Ka emission 
line, setting cr=0 to simplify the fit. The best fitting energies and equivalent widths for the emission line 
are shown in Table |2l The line energy is consistent for all the datasets, although the intensity of the line 
does change with the AGN state and the local absorption column, as expected. The errors o n the equivalent 



widths are quite big, which may be a sign of the line having some breadth, as suggested by iFeldmeier et al 
(|l999h . or be a side effect of the complicated model. 



To compare our results with the previous ASCA and XMM results, and verify our hypothesis of variable 
absorption, we tested a few other models. We found that a single p ower law, with local absorption, or local 
+ intrinsic absorption, results in values similar to those obtained by llmmler et all (l2003h . although it gives a 
very poor fit to the 2003 XMM spectra and does not fit at all the Chandra data (reduced ~6). This is most 
likely caused by the varying absorption column, which cannot be properly accounted for with such a simple 
model. 

We found that the only model that could allow a direct comparison with all of the previous results 
was one consisting of local absorption, a double partial covering (zpcfabs), a power law and a Gaussian 
to account for the Fe emission line. We tied together the covering fractions across the different spectra, 
allowing only the Nh to vary. We did not add any reflection components to this model, since they are 
not considered in the ASCA analysis by IFeldmeier et all (| 19991) . but we did add the contribution due to the 
thermal emission from the radio bubbles. The results of this fit are shown in Table|2]under the label "Model 
2". We find that this model gives 3046 for 2978 degrees of freedom, reduced ;^'^=l. 023. We found that 
the photon index of the power law i n our fit is lowe r than the one obtained in previous analysis of Mrk 6, 
only similar to the value obtained by IS church et all (|2006') with their third model. This is partly due to the 
addition of the thermal emission component, which we deem necessary since it has been spatially resolved 
and fitted to the Chandra data and results in a statistical improvement in our other model. It is also partly 
caused by the effect of the variable absorption. We find the values of the Fe Ka emission line energy and 
equivalent width to be consistent with the ones we obtained for the previous model. We found that it is 
necessary for the two partial covering components to overlap ( f!^„„ + f}„„ > 1) for th e fit to be statistically 
acceptable. We tried to fix the values to the ones employed by ISchurch et al.l (|2006h in their analysis, but 
this results in a reduced ~1.7. 



While llmmler et all (120031) cite a value for this photon index of 1.81 from their analysis of the Beppo- 
SAX data, which is less affected by variable absorption, they mention in their analysis that t he source was 
in a hi gher state when those data w ere obtained. Anothe r Beppo-SAX observatio n analysed by lMalizia et al. 
(l2003h gives a photon index -1.5. ISchurch et al.l (120061) : IFeldmeier et al.l (Il999l) also favour a higher value 
for the photon index (1.6-2.2), but setting the photon index to 1.6 results in a reduced ~1.13 in our joint 
fit, and the fit requires a heavier absorbing column for the second partial covering (Nh = 66.61^4 ^4 x 10^^ 
cm~^ for the Chandra spectrum). 



As shown in Table [2l the absorption column is over an order of magnitude laiger in t he Chandra data 
than in the XMM, for both models. This high value is more similar to the ones obtained by IFeldmeier et al. 
(I1999I) in their analysis of the ASCA data. The hypothesis of "obscuring" clouds near the AGN was proposed 
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by llmmler et alj (120031) . and lSchurch et all (120061) supported it after obtaining a relatively smaller absorption 
column in their analysis of the XMM data from 2003. This detection by Chandra rules out instrumental bias. 
Although more data are needed to make a statistical test, current evidence seems to point towards denser 
material having been observed in 1997 (ASCA) and 2009 (Chandra), with the passing clouds being least 
dense in 2003 and intermediate obscuration values in 2001 and 2005 XMM-Newton. This allows us to 
establish a variability timescale of 2-6 years, which translates to a distance to the central black hole of 
just a few pc. Since grating data would be very useful to get an additional insight on the structure and 
composition of this variable absorption, we briefly investigated the XMM-Newton RGS spectra from the 
three observations. Unfortunately, the signal to noi se is rather bad, and no e vident emission or ab sorption 
lines can be inferred from them, as already noted by lSchurch et al.l (120060 and llmmler et all (120030 . 



Although the statistical result is better for the second model, we must consider the physical implications 
of both models. This AGN seems to fall between the behaviour expected from a classical Seyfert and that 
of a radio galaxy, it is therefore interesting to explore the consequences of th ese results in the light of what 
we know about these classifications. As explained in lHardcastle et all (|2009l) . the soft excess characteristic 
of most Narrow-Line Radio Galaxies (NLRGs) can be explained by three classes of model: 

• It is thermal or line emission, from the photoionized material close to the AGN or from the IGM of 
the host galaxy. 

• It is non-thermal, power law emission from the central AGN, visible through the partial covering 
material. 

• It is non-thermal, power law emission related to the jet. 



The first model can most likely be discarded in this case, since neither the Chandra nor the XMM 
spectra show evidence for strong emission lines. The RGS spectra do show some residuals over a pur e 
continuum model in the 0.5-1 keV region that could be a hint of oxygen emission lines lSchurch et al.l (120061) . 
but they are quite faint, and since the area is not spatially resolved, they could be due to the thermal emission 
from the radio lobes. 

To test which of the two other models is most likely to apply in this case we calculated the X-ray 
luminosity of the unabsorbed power law from our first model, Lx^u = 9.6 x 10^^ erg s"^, and the vFy 
5G Hz luminosity at the ba se of the jet from our radio maps, L^ghz. - 1-41 x 10^** erg s~^ On the plots 



by iHardcastle et al.l (l2009h this falls quite far from the behaviour expected for radio galaxies. We must 
therefore conclude that the jet-related soft emission must be quite low, and the dominant contribution to the 
soft excess is most likely the emission from the central AGN. 



Following the results of Evans et al.l (120101) we also applied their model to our data, to check if it 
would fit a photon index for the power law more consistent with the values expected from other inter- 
mediate Seyferts. This model tests a different absorption component for the soft and hard power laws: 
wabs(apec+zwabs*po+zwabs*z.pcfabs(zgauss+po)). We obtained a fit to the data similar to that of our pre- 
vious partial covering model, =3037.22 for 2971 DOF, reduced ;^^= 1.022, and a similar photon index for 
the second power law, 1.25. 
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We also tested a model with a multiphase medium. We found that a single warm absorber does not 
provide a good fit to the data, and hence split the intrinsic partial covering absorptions of our second model 
into a cold (zpcfabs) and a warm (zxipcf) component. This model seems to provide a better fit for the soft 
end of the spectra, and yields a;^'^=3000.98 for 2977 DOF, reduced ;^'^= 1.007. The absorbing columns are 
similar to those we found for the double cold partial covering, and the photon index of the power law is 
1 -41^003, higher than with any other model we have tested. However, the ionization parameter is not very 
well constrained, logE = -O.OQ^JH!, since the model is limited by the resolution of ou r spectra. Warm 



absorbers have been observed in other Seyferts, albeit at higher ionization values (see e.g. iLonginotti et al. 



20091) . but in our case the situation is less clear. The soft excess this model successfully accounts for can 



have several origins: the jet, warm-hot gas from other regions of the galaxy, or a higher normalization of the 
thermal emission from the lobes, since we decided to freeze the latter to the maximum value provided by 
the spatially resolved fit. Unfortunately the grating spectra are not sensitive enough to test either hypothesis, 
but future, longer exposures may hold the key to the complexity we are observing in Mrk 6. 



3.2. Radio bubble-related emission 

To extract the spectra of the extended radio bubbles we only used the Chandra data, since in the 
XMM images this area is com pletely dominated b y the PSF of the central AGN. We used the contours from 



the radio images obtained by Kharb et al.l (120061) to define our extraction regions around the 7-kpc radio 
bubbles, which can be seen in Figure [T] The photon statistics for the bubbles are not very good; we do not 
have enough counts to fit to both bubbles separately and determine their temperatures independently, and 
even for the joint extraction the photon count is quite low (~200 counts after background subtraction, see 
Figure IT). Even though we were careful to keep these regions well away from the central AGN, there is 
still some contamination from its PSF in our spectra. The data we simulated with ChaRT and MARX to 
correct the effects of pileup in the AGN spectrum (see Section [3?T]) show the number and energy distribution 
of photons from the AGN PSF across all the CCD, see Figure [3l and can thus be used to estimate the AGN 
photon contribution to our spectra from these external regions. Our calculations show that this contamination 
amounts to up to 60% of the counts in the spectra prior to background subtraction. However, the background 
regions that we have chosen contain a similar number and energy distribution of photons from the AGN, its 
PSF being quite homogeneous on these scales, and their effect is hence mostly eliminated. We carried out 
a study of the AGN contamination in our regions by studying the energy distribution of the AGN counts in 
our simulated data, within the same extraction regions. We found most of the counts to be contributing at 
energies above 4 keV. There is a minor contribution at energies slightly above 1 keV. The residuals we get 
from our fit at these energy ranges are likely to have an origin in this contamination (see Figure IT). 

We have found a source of contamination from high energy photons in the form of a small bright region 
(Lx ~ 1.5 X 10^^ erg s~') that can be seen on the left panels of both Figure[T]and Figure[3]just South of the 
source (labeled as "A" in the left panel of figure [Jl. We plotted a histogram of the counts from this region, 
and by their distribution, and the fact that it has no radio or optical counterpart, we infer this is clearly a 
point source, most likely a background, more distant quasar. It is quite unlikely that this object may be an 
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X-ray binary (XRB), but it could be an ultra-luminous X/ray source (ULX), since its magnitude should be 
above our detectability threshold for this distan ce {Lur.x ~ 10^^ - 10^^ erg s^^) although the morphology of 
the host galaxy indicates this is not very likely dSwartz et alj|2004l) . 

The other bright structure North of Mrk 6, labeled as "B" in the left panel of Figure [H does seem to 
have a thermal origin, and a similar energy histogram as the photons we find in the bubbles; part of it is 
clearly surrounding the radio structure, and we therefore assume that it is part of the hot shell outside the 
radio bubbles and include it in our extraction regions. It must be noted, however, tha t Mrk 6 has an inn er set 
of ~1.5 kpc radio bubbles, almost perpendicular to the outer ~7.5 kpc structure (see Kharb et al. 20061) . and 
this is likely to cause some turmoil in the gas surround ing the AGN. These inner bubbles are not resolved by 
Chandra. The jet of Mrk 6 was detected by MERLIN (IKukula et al.ll 19961) and is aligned in the North-South 
direction, so that this is the presumed photoionizing axis of the AGN. 

This direction also c orresponds to the Ex t ended Narrow-Line Region (ENLR) of Mrk 6, see Figure |2] 
and ' Capetti et all d 1995b : kukula et aP (Il996h : llGiarb et all hood) for details. An overlap between [OIII] 
and X-ray emission has been observed in many Seyfert 2 AGNs before, and it is believed to happen when 
the X-ray photons heat up and ionize the cold gas around the AGN. The [OIII] emission in this system 
seems to correspond to different structures: the emission closer to the active nucleus is clearly caused by 
direct photoionization from the AGN, while the two bulges that extend in the N-S direction, showing both 
enhanced [OIII] and X-ray emission, are most likely caused by the interaction of the current jet with the 



surr ounding material, which heats it up both through shocks and photoionization (see ^ Capetti et all 1 1995 



and Kukula et al.ll 19961 for details on these inner structures). Since the overlap between the [OIII] and X-ray 
emission can also be observed at kpc distances, extending to the North of the AGN, we calculated whether 
the nucl eus could be photo ionizing these regions to X-ray and [OIII] emittin g temperatures, following the 
steps of IWang et all (120091) and the models from iKalhnan & McCrayl (119821) . For the X-ray emission that 
can be seen ~3 kpc North of the AGN, around the edge of the radio structure, we used the parameters 
derived from our fits of the nuclear spectrum and the radio lobes. We found that such X-ray emission cannot 
be caused by photoionization from the AGN, since the unabsorbed luminosity of Mrk 6 (~ 1.3 x lO'*-' erg 
s"') can only account for photoionized X-ray emission at distances within 300 pc of the central source. The 
[OIII], however, could be photoionized by the AGN at much larger distances, as long as the electron density 



stays low (< 10^ cm^ for distances over 3 kpc). It is likely that some of the [OIII] emission to the North 
of the source may have been caused by the shock itself, both via direct heating and photoionization. The 
fact that the [OIII] emission follows the outline of the radio structure, but is mostly detached from it, also 
supports this scenario. 

With our photon statistics we cannot rule out a model in which this emission is produced in some way 
by a precessing jet, which would also explain some of the emission in the ENLR. Both assumptions could 
explain the structures that surround the real, but not the simulated AGN, on Figu re |3l For de t ails on the 
radio structure of Mrk 6 and the discussion of a possible jet precession scenario see lKharb et all (|2006r) . 

We fitted to our data an emission spectrum from collisionally-ionized diffuse gas (apec model in 
XSPEC) with the abundance fixed at 0.35© and Galactic absorption, using again the averaged column den- 
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Fig. 7. — Spectral fit to the spatially constrained radio bubble regions, with the best fitting apec (T=0.94_Jj jg 
keV, abundance=O.350, ;i-2=9.39 for 8 PHA bins and 6 DOF). 



sity result Nh - 6.39 x 10^" atoms cm"^ from lDickey & LockmanI (Il99d) . We estimated a temperature for 
the bubbles of 0.94^^ }^ keV. This model yields nx^ value of 9.39 for 8 PHA bins and 6 degrees of freedom. 
Different extraction regions yield slightly lower values of kT , but still well within the errors, and different 
values oix^, due to the reduced number of bins. \ ye also tested a l arger abundance of 0.6© consistent with 
our estimations on the thermal emission of Cen A dlCraft et allboosl) . but this requires a lower normalization 
and yields a poorer fit to the data (^^ of 11.42 for 8 PHA bins and 6 DOF). 

The normalization of this component is 8.5;^^} ^ x 10~^cm~^, just high enough to be noticeable on the 
spectra we extracted for the AGN (~2% of the flux between 0.8-1.1 keV). We added an apec component to 
our AGN model and performed a statistical test on the temperature, finding it to agree with the results from 
the spatially constrained spectrum (^r = 0.87;';g 23 keV). This addition also improves the statistics of the 
joint fit, if only slightly (reduced reduced from 1.051 to 1.049). The effects of adding this component to 
the Chandra spectrum can be seen on Figure |5] 

To discard the possibility of the emission being non-thermal synchrotron or inverse Compton, we at- 
tempted to fit a model consisting of a power law and local absorption to the data, and obtained very poor 
results (reduced of 3.79). Therefore we can quite safely assume the emission to be thermal. 

Although the abundance that best describes the thermal emission from the bubbles is 0.35©, more 
consistent with NGC 3801 (0.3-0.4©) than with Cen A (0.6©), the statistical difference alone is not enough 
to rule out a higher value for the abundance, due to the low number of photon counts. We could argue 
that we are studying a quite small early-type SO galaxy which probably has, as our results show, a very 
low gas density outside the central regions. The dust lane that is hinted at in the optical images may be 
a result of a past merger, there is an equivalent structure in NGC 3801. Furthermore, we must consider 
the influence of the active nucleus on the star formation liistory of Mrk 6. While it has been argued for 
more po werful systems that the AGN activity may indeed be s uppressing star formation in the host galaxy 
(see e.g. iMcNamara & NulsenI I2007L ISchawinski et al.l l2009l iNesvadba et all 12010 ) the situation may be 



-17- 



more complex for lower power systems, and early type galaxies in particular, as suggested by ISchawinski 
(I2010I) . While these arguments might be enough to support the low abundance we are observing, our tightest 
constraint comes from imposing the shock conditions for the bubbles. The details of these constraints and 
their consequences will be discussed in Section 



Even though the data show a significant correspondence between the radio contours and the X-ray 
emission inside the North bubble (see Figure[T|l, which could be associated with jet/ISM interac tions leading 
to the heating of gas in these regions, as in NGC 6764 dCroston et al.ll2008al iKharb et al.ll2010h . most of the 
X-ray photons seem to correspond to the edges of the bubbles or even the regions just outside them, which 
is mo re compatible with t he presence of shells of strongly shocked gas outside the bubbles, just as in NGC 
3801 dCroston et al.ll2007h . While some structure can be seen both in the radio maps and in the X-ray image, 
we will assume in a first approach that the shells are spherical and uniformly filled with gas, and will later 
discuss the implications of a non uniform density. 

The fact that the radio emission is apparently fainter in the South bubble (see Figure [B, and the small 
difference in the apparent sizes of the bubbles, the North one being slightly bigger (see the discussion of the 
size of the shocked shells in Section l34l ). is probably due to a different past history for both bubbles, due to 
the gas expanding in each direction having encountered a surrounding medium with a different density, or 
to an asymmetric energ y output from the AGN. The radio images show there is also a similar asymmetry in 
the inner E-W bubbles ( Kharb et allbooeh and the N-S current jet j xukula et aDll996[) . The [OIII] emission 
is also asymmetric (see Figure |2ll. These asymmetries are not visible in the optical (HST) images of the host 
galaxy, the apparent dust lane being roughly perpendicular to the outer radio bubbles, but if the jet is indeed 
precessing, it could have created a somewhat cluttered environment around the AGN. A precessing jet is 
not likely to explain an asymmetry in the energy fed by the AGN to each bubble, but we cannot rule this 



out. D oppler boosting has been suggested for some Seyfert jets such as the one in Mrk 23 1 ([Reynolds et al. 



20091), but it is unlikely to be the cause of the asymmetry of the lobes in Mrk 6, since it would require 
relativistic speeds for the gas in the shells, which is inconsistent with our results (see Section 13.41 ). This 
asymmetry only underlines how challenging it is to make a self-consistent model for this kind of object. We 
will assume that these effects are minor for our purposes. 



3.3. Luminosity constraints on the external environment 



To establish whether the gas from the bubbles is driving a shock into its surroundings we need to 
estimate the temperature of the external r tiedium, derive the relati ve pressures and impose the Rankine- 
Hugoniot conditions for pressure balance (ILandau & Lifshitzlll987n . However, we have found that due to 
the galaxy's morphology and size, the halo regions where the gas is expanding are not bright enough in the 
X-rays to allow us to obtain a spectrum. While the Chandra spectrum appears to admit an additional thermal 
component with kT ~0.35 keV, there is no statistical improvement on the XMM spectra from adding it; we 
therefore cannot state that the component is real. 



To obtain an upper limit on the luminosity and normalization of the gas component in these halo 
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Fig. 8. — Lx vs. Lb plot from 10' Sullivan et all (120011) with superimposed Mrk 6 B band bolometric lu- 
minosity (vertical blue line) and the expected (upper limit) X-ray luminosity range we calculated for the 
extended thermal emission outside the shells (horizontal green strip). The diagonal dashed line is the best 
fit to the data excluding AGNs, BCGs and dwarf galaxies. The AGN is most likely causing a blue excess in 
the optical, so the B-band luminosity is an upper limit. Given that the temperature of the external medium 
is most consistent with 0. 1-0.2 keV, an X-ray excess caused by this emission is likely to be happening. 



regions, we calculated the background-subtracted counts from a wide annular shell centred in the AGN 
{Rint=2l arcsec, Rext=59 arcsec, excluding a few point sources resolved by Chandra). We calculated the 
relative contribution of counts from the AGN using the simulated dataset, applying a scaling factor obtained 
from the ratio of counts within the innermost regions of the PSF (excluding the pixels we used for pileup 
correction) in the real and simulated data, which takes into account possible energy-dependent variations. 
We then subtracted the AGN contribution from the total number of counts in the halo, and added the 3cr 
error, obtaining an upper limit of 270 counts, equivalent to a rate of 3.65x10"^ counts s"^ 



-19- 



Using these statistical constraints on the number of photons from the halo, we can derive upper limits 
on the luminosity of this region. These upper limits also constrain the possible temperature, electron density 
and pressure for the gas, allowing us to test the shock conditions. To achieve this, we chose a range of 
possible temperatures for the medium outside the radio bubbles (0.1-0.6 keV) and used XSPEC with a toy 
apec model to determine the limit on the normalization corresponding to our count rate for each of our 
chosen temperatures. We maintained the Nh and redshift parameters fixed to the appropriate values for the 
system (see previous sections), and the abundance fixed at 0.35©; we will discuss the implications of larger 
abundances in the next Section. We then derived from the flux in these models the range of possible X-ray 
bolometric luminosities (see Table [3]). 



We then estimated the B band luminosity of Mrk 6 from the data in the catalogue by lde Vaucouleurs et al 
(11995 b. obtaining a value of LISxIO^^LbO. Since Mrk 6 contains an AGN, it is likely that we are observ- 



in g a blue excess, therefor e this value is to be considered an upper limit. We used as a reference the work 



of 10 'Sullivan et all (120011) to estimate the feasibility of our results on the X-ray luminosity. Their work 
catalogues the X-ray extended emission detected by ROSAT on several hundred early-type galaxies, and at- 
tempts to establish a relationship between the X-ray and B band bolometric luminosities for these systems. 
We overplotted our results on their diagranj^ (see Figure [8]l to see whether Mrk 6 falls near or above their 
correlation. They excluded AGNs on their fit, and are thus only measuring the extended emission from the 
galaxy, which is what we are looking for in the outer regions of Mrk 6. Our results cover a wide range of 
possible X-ray luminosities, but especially for the lower temperatures, and considering that our estimate for 
the B luminosity is an upper limit, the data do fall in the right region of the plot. We can therefore assume 
that, although we are not directly resolving it in our Chandra image, our upper limit of counts is consistent 
with the luminosity expected for the galaxy, excluding the AGN. This means these photons indeed come 
from the extended X-ray emission of the galaxy, not from the background or the AGN. 

We also tested the K-band luminosity of M rk 6 (l.lTxlO^^L©) against its X-ray luminosity using the re- 
lations proposed by Mulchaey & Jelternaf (1201 0). to avoid the effects of blue excess fro m the AGN. We found 



the res ults to agree with our estimations from the B band. Our results also agree with iMulchaey & Jeltema 



(|2010h on the hot gas content of a sample of early-type galaxies in different environments. Their work shows 
that although some of these objects seem to retain their hot gas halos when living in clusters, others lose it, 
while for isolated galaxies the hot gas content is mostly related to the mass of the galaxy. Mrk 6 is quite 
small and seems to be living in a poor environment; it is therefore likely to have lost most of its halo gas. 

Finally, it must be noted that in the ChaRT/MARX simulations the wings of the PSF are known to 
be underestimated at large distances from the point source, especially at low energies. We contacted the 
Chandra Calibration team via the HelpDeslH and consulted the available bibliograph>0, and estimated that, 
for the regions we are considering in the halo of Mrk 6, the simulated PSF wings could be underestimated 



^We obtained the data from |http://www.sr.bham.ac.uk/~ejos/catalogue.htnil| 
' http ://cxc .harvard.edu/hef pdesk/ 

''See |http://cxc.harvard.edu/cal/Hnna/hrma/psf/index.htmf I 



by up to a factor of 2. If this were the case, our upper Hmit of counts would be more conservative (~30% 
lower) and we would be near the limit of non detection, implying an external gas temperature below 0.1 
keV (see discussion in the next Section). If the gas were much colder than 0.1 keV we would not be able to 
constrain any limits on the density or assume a shock scenario. This would not be very c onsistent with wha t 
we know about the virial temperatures of the halo gas in early type galaxies (see e.g. iNaab et al.l 120071) . 
but since we do not have the instruments to measure the extreme UV emission a gas of these characteristics 
would produce, it cannot be ruled out. We must also note that the X-ray emission cannot be accounted for by 
photoionization mechanisms (see discussion on the [OIII] emission in the previous Section). We conclude 
that the shock scenario provides the most plausible explanation for the emission we are observing. 



3.4. Physical properties of the gas shells and ISM 



Table 3: Estimated apec normalization, upper luminosity limits, electron densities and gas pressure outside 
the bubbles for a sample of possible gas temperatures. 
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Table 4: Top: Model normalization, estimated upper and lower limits for the shells, derived limits for rie^out required to comply with the 
Rankine-Hugoniot conditions for a strong shock, and inferred limits on the mass, pressure, total thermal energy, wo rk available from th e 
gas filling the shells and total kinetic energy. Bottom: Same estimations for a set of Cen A-like 200 pc-wide shells dCroston et alj|2009l) . 
The errors on the pressures and energies are derived from the errors in the temperatures of the shells in the XSPEC fit (kT = 0.94+°- }2 keV) 
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Using the results from the previous section we can now calculate the upper limits on the electron 
densities, n^, for each temperature, from the definition of the model's normalization: 

10-14 

fieNndV (1) 



14 n 
+ Z)f J 



47r[DA(l 

where Da is the angular diameter distance to the source. We assumed the volume to be the one mentioned 
in the previous subsection, a wide spherical shell centred on the AGN and beginning just outside the edge- 
brightened emission around the radio bubbles. The results are displayed in Table [3] Notice how the derived 
ne,out increases sharply below temperatures of 0.3 keV. This is caused by the loss of sensitivity of the instru- 
ment at low energies, and a change in the behaviour and type of emission lines found at these temperatures. 
We can derive from these values the pressure of the gas, also displayed in Table [3] The pressures are con- 
sistent with the environment of Mrk 6 being quite cold and rarefied. Notice the increase of the pressure for 
very low temperatures: this is a consequence of the steep increase in the electron density estimations. Our 
choice of the ~60 arcsec region is based on the most likely dimensions of the halo; in order to test these 
constraints we selected a larger region arcsec) and found that the densities and pressures change 

by 15% at most. 

We also estimated the limits on the electron densities for the shells, to compare them with these upper 
limits we obtained for the surrounding gas. This is necessary to verify that the proposed shock scenario 
holds, and see whether any limits on the external temperature can be derived. To do this we assumed that the 
dominant thermal contribution to the spectral fit we obtained for the radio bubbles (Section 13.21 ) originates 
in the edge-brightened, shell-like emission around them. Our extraction region covers the whole surface of 
the bubbles to include any photons from the nearest and farthest areas of the spherical shells, which have a 
lower surface brightness due to projection effects. The gas inside the bubbles is radio emitting plasma, and 
therefore does not contribute to our spectrum, and any emission from gas immediately outside the shells is 
accounted for by the background. 

Using the constraints on the temperature and model normalization from our bubble fit, we then es- 
timated the North and South shell thickness from the apparent thickness of the edge-brightened emission 
around the radio bubbles, which in Mrk 6 is quite evident, see Figured] Since the spectral resolution is not 
enough to allow us to discern the detailed structure of the shells, we assumed a simple spherical geometry 
for each structure. From the extent of the emission we estimated for the North shell a width of ~980 pc 
{Rint,N = 2.24 kpc, Rext,N = 3.22 kpc), and for the South shell -1100 pc {Rint,s - 1-55 kpc, Rext,s = 2.65 
kpc). Estimating the thickness of the shells is difficult due to the poor photon statistics, hence these values 
were taken as a conservative upper limit. They are larger than what we expected a priori, but this may partly 
be due to the apparent lower density of the external environment in Mrk 6, as well as to the limited reso- 
lution we can achieve with our photon statistics. The results we obtained are shown in Table |4] This table 
also shows what the limits on the would be, both for the shells and the external environment, assuming a 
thickness of ~200 pc, similar to the one estimated for Cen A. 

From these values we can derive the total mass of the gas contained in each structure and the resulting 
pressure. These res ults are also displayed in Table |4] When comparing these results with NGC 3801 
( Croston et al we find that even though the electron density is lower in Mrk 6, the apparent thickness 
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of the shells is much bigger, hence the total mass of the gas is higher by slightly less than one order of 
magnitude. The inferred pressures for the gas are roughly the same as the ones found on NGC 3801, due 
to the temperature of the gas being higher in Mrk 6. The errors on the pressure are calculated from the 
uncertainty in the gas temperature, 0.94^°- keV. We can see that choosing thinner shells, like the ones 
found on Cen A, would imply a much smaller (~0.5) fraction of gas, while the pressure would almost 
double (see lower section of Table H)). The pressure jump Psheii/Pout for a 1 kpc shell is ~7-63, higher than 
in NGC 3801, but again this is to be expected, due to the bigger temperature contrast and more rarefied 
external environment in Mrk 6. The contrast would be much higher for thinner ~200 pc shells, ~ 14- 122. 

We calculated the minimum internal pressure of the Northern bubble under equipartition conditions 
from the radio data, fitting a broken power-law electron energy spectrum with p = 2 at low energies, 
steepening to = 3 at the electron energy that gives the best fit to the data. We ass ume no protons an d 
7,„m-10. We obtained a value of P~4xl0~'^ Pa, consistent with the results obtained by lKharb et all (|2006r) . 
and roughly an order of magnitude lower than the pressures we deriv e from our X-ray dat a. This depar- 
ture from minimum energy is often found in FR I radio galaxies (e.g. iMorganti et al.lll988h . and although 
there are very few examples where the minimum internal pressure is higher than the external pressure in the 
ISM, we also found this effect in NGC 3801 and Cen A. These results imply that there is some additional 
contribution to the internal pressure, caused by a large deviation from equipartition conditions, a high frac- 
tion of non-radiating particles, such as thermal or relativistic protons originated from the interaction of the 
galaxy's gas with the j et, or a low filling fac tor. Of these, the second explanation is the most plausible (see 



Hardcastle et al 



Kharb et al 



20071 ICroston et al.ll2008ah . even more so if the jet is indeed processing, as suggested by 



(|2006r) . or if it is just episodic. 



In the limit where there is no shock t he pressures of the external medium, P„„,, and the shells, P shell 
should be equal (ILandau & Lifshitzl 119870 . The difference in the apparent luminosities of both regions 
already rules out this situation if both temperatures are equal. We can also rule it out by noticing that these 
values of «e,„„; and P„„, would imply external temperatures far below 0.1 keV (see the results from Tables [3] 
and HI) which is unlikely, as discussed in the previous Section. The density ratio Poutip shell tends to 4 in the 
case of a strong shock. For values of poutIp shell smaller than 4 the required values of the external temperature 
are below or just above 0. 1 keV; the most plausible scenario is therefore that of a strong shock. 

Assuming a density contrast of 4, and using again the limits on Ue^out from Table[3l we can constrain the 
possible temperature for the external gas to be around 0.2 keV. This result is consistent for both the North 
and South shell values and for both the 1 kpc and 200 pc shell thicknesses (the 200 pc hypothesis would 
imply a higher gas density within the shells, thus to maintain the density contrast of 4 the upper cap on the 
external temperature would be lowered, but the results still imply kT > 0.\ keV). We were able to directly 
measure the temperature of the external m edium in Cen A (IKraft et al.ll2003f) . and found the temperature of 
the ISM to be -0.29 keV. For NGC 3801 Jcroston et al.ll2007h we obtained a value of 0.231°-^^ keV. These 
values directly applied to Mrk 6 would cause an inconsistency, since the density ratio cannot exceed a value 
of 4, but given the uncertainties they are still in very good agreement with the current results. 



As an additional test we repeated our calculations for Z=0.6o and Z^l.Oo. The model normalizations. 
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and consequently the electron densities, inversely depend on the abundance, so we found lower values as 
we approached the solar values. However, the impact of this effect is bigger on the results from the external 
gas than to that contained in the shells. This causes an inconsistency when applying the conditions for a 
strong shock, since the external gas would be required to have temperatures far below 0. 1 keV, which are 
implausible, as we discussed above. Higher abundances also have an impact on the energy yield. The 
pressure and total energy for Z=0.6o would be ~ 20% lower than the ones we calculated with Z=0.35o, and 
~ 40% lower for Z=1.0o. 

To further test these res ults, we used the B-band lurn inosity of Mrk 6 to relate it to other galaxies 
following the relations from lO'Sulhvan et al.l (120011 . |2003|). We found that both the temperature of the 
external gas (0.2 keV) and the abundance we fit to the data (0.35©) are consistent with their relations. Their 
sample includes two galaxies with gas X-ray and B-band luminosities slightly larger (Lx ~ 10'*° erg s~^ 
Lb ~ IO'^Lq) than what we found for Mrk 6 where these parameters can be directly related. For NGC 1549 
they find an outer gas temperature ~0.25 keV and an abundance -0.140. For NGC 4697 their results show 
T ~0.24 keV and Z ~0.4o. We therefore consider our best estimates of the properties of the large-scale hot 
gas in Mrk 6 to be plausible ones. 



We can calculate the Mach number for the shock, as follows (ILandau & Lifshital 19871) 



4(r+i)(r,/„,,/r,»,) + (r-i) 

M = y (2) 

where F is the polytropic index of the gas, assumed to be 5/3. The results for our chosen temperature range 
are displayed in Table [S] 



Table 5: Possible values of the Mach number given by our temperature constraints. The values of T shell are 
the best fit value and the IcrUmits allowed by the errors. The most likely value is shown in bold. 



kT„,„ (keV) 


kTsi„ii (keV) 
0.75 0.94 1.06 


0.10 


4.92 


5.50 


5.84 


0.20 


3.49 


3.90 


4.14 


0.30 


2.86 


3.20 


3.39 


0.40 


2.49 


2.78 


2.95 


0.50 


2.24 


2.49 


2.64 


0.60 


2.05 


2.28 


2.42 



From the previous constraints we can assume the most likely value for the temperature to be 0.2^q j 
keV. This yields a value for the Mach numbe r of 3.90||;} q4, of the same order of magnitude as the one we 



obtained for NGC 3801 (ICroston et al.ll2007r) . although smaller than the one found for Cen A dKraft et al. 



20031). This is to be expected, since the morphology of Mrk 6 is more similar to the one of NGC 3801 than 



to that of Cen A, where only one shell is detected. 

The speed of sound in the ISM ranges from -190 km s~' in the case where the external temperature 
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is 0.1 keV to ~460 km when kT = 0.6 keV. Since the Mach number also varies with the temperature 
(see Table [5]), we derived a consistent value of 1030;|;jjq km s"^ for the speed of the gas. This allows us to 
derive the kinetic energy of the gas, which is displayed in Table \4\ The total thermal energy and total work 
available for each radio bubble are also displayed in Table |4l 

The total energy (thermal + kinetic) from both radio bubbles, under the 1 kpc width assumption, is 
~ 2.6-4.6 X 10^^ erg. This corresponds to ~3.6xlO^ superno va explosions with individual explosive energy 



of 10^^ erg. This result also agrees with the value derived by iKharb et all (I2006D from the radio data of Mrk 
6, ~ 1 X 10^^ erg, since it was a lower limit. They estimate a kinetic luminosity for the jet of ~ 10^^ erg s~^, 
derived from an inferred radio luminosity of ~ lO^'' erg s~' and an efficiency of 1%. With our results, the jet 
would require ~ 1.1 x 10^ years to deposit all its energy into the nuclear ISM and produce this pair of radio 
bubbles. This result is almost two orders of magnitude larger than the one they obtained, but not incongruent, 
since their estimation was a lower limit derived from the equipartition pressure. These results also agree with 
the timescales they propose in the context of a precessing jet. Another timescale can be inferred from the 
supersonic expansion speed of the bubbles, which would require a time of~2.5-3.1xlO^ years to inflate 
to their current size, assuming M ~ 4 throughout their lifetime. If we adopt this timescale, the required jet 
power is ~ 7 X 10"^^ erg s~^ 

The thermal energy stored in the shells is an order of magnitude larger than the kinetic energy required 
to inflate the cavities (see Table IDl, meaning shock heating must be the dominant energy transfer mechanism 
during this stage of the evolution of the radio source. Although the emission from the surrounding gas is 
too weak for us to fit a model to its brightness profile, which would be ideal to determine the direct impact 
of this thermal energy into the surrounding ISM, our results, and comparison with those obtained for Cen 
A and NGC 3801, show that this energy is enough to significantly alter the surrounding ISM, sweeping out 
the remaining halo gas. 

With our current photon statistics we have no direct evidence of changes in the brightness (and therefore 
gas density) across the shells, therefore we cannot state whether these structures are brighter in the inner 
edge, closer to the nucleus, in Mrk 6. However, this would be expected if the bubble expansion is supersonic 
everywhere: even though the shock is weakest close to the nucleus, the ISM is much denser in this direction, 
hence the associated X-ray emission is higher. 

If we assumed the density of the gas not to be constant across the shell, with a higher density in the 
area closer to the radio source, the derived constraints on the temperature of the external medium would be 
relaxed, allowing higher values. Consequently, the pressure and density of the external gas could also be 
higher. This situation would bring our results closer to the values we obtained for Cen A and NGC 3801. 
While this would be desirable for consistency, we must note the scenario in Mrk 6 is quite diff"erent, both due 
to its morphology and to its lower radio power, thus we could indeed be seeing a case of the AGN driving 
a strong shock into a very cold environment. It could also be the case that the shells are not continuous all 
around the radio bubbles, which would change the constraints on the density of the gas. This hypothesis 
could be tested with a longer exposure. 

To conclude our analysis and study the fueling mechanism of the AGN, we used our estimates on the 



external gas temp erature and density to ca lculate the Bondi accretion rate from the central hot gas. We used 
the relations from lMarconi & Hung (120031) and the K-band luminosity of Mrk 6 and obtained a mass for the 
black hole of S.VxIO^Mq. The Bondi rate is given by 



M = AnApG^Ml^cf (3) 

where A=0.25, p is the density of the gas at the Bondi radius, and cs is the speed of sound in the medium 
(cs = ^jTkT //^nip). Since the density and temperature at the Bondi radius cannot be lower than the ones we 
estimated for the gas outside the shells, we used those values to obtain a lower limit of M = 2.04 x lO'^M© 
yr"^. Assuming ju=0.l this translates to a Bondi power Pb = 1.1 x 10"^^ erg s"'. The energy output we 
measured from the AGN bolometric luminosity (~ 2.4 x 10^^ erg s"^, using a bolometric correction factor 
of 19) and the shocked gas shells (~ 7 x 10^^ erg s~^) imply an accretion rate of M ~ 5.7 x 10"^ Mq yr"^. 
This suggests that Bondi accretion is unlikely to provide all the energy we observed to be transferred to the 
ISM, even with an efficiency of 100%. Assuming /i=0.1, the density would have to be over three orders 
of magnitude higher than our estimates of the density of the unshocked gas in order for Bondi accretion 
to be viable. If Bondi accretion cannot account for the energy we are observing, the main fuel source for 
this system is likely to be cold gas. Cold accretion does not depend on the black hole mass and can happen 
in poor environments as long as there is a fuel reservoir, which may be the case in Mrk 6 if the central 
obscuring material seen in the optical is indeed cold gas from a past merger. 



4. Conclusions 

We conclude that there is a high probability that we are observing shells of strongly shocked gas around 
both radio bubbles of Mrk 6, with Mach numbers 3.2-5.5 consistent with the Rankine-Hugoniot conditions 
for a strong shock. This is the third clear detection of such a process in a low-power radio source, and 
therefore reinforces the hypothesis that this mechanism may be very common and may play a fundamental 
role in the process by which these young radio sources form their stellar populations. 

It is very likely that much of the gas in the shells will escape the gravitational potential of the host 
galaxy. Our results show shock heating is the dominant process, thermal energy being an order of magnitude 
larger than the work required to inflate the radio cavities, implying the outburst impact on the host galaxy's 
ISM is likely to be dramatic. Moreover, the inner set of bubbles is expanding perpendicularly to the external 
ones, and into denser regions of the host galaxy. Although the shock will probably be weaker, due to the 
increased density, the expansion of the inner bubbles is likely to cause a more direct disruption on the 
ISM and are likely to have a bigger impact on star formation. Unfortunately the resolution of the current 
generation of X-ray telescopes does not allow us to determine the details of this ejection of gas. 

We also conclude that the hypothesis of a variable absorption gas column, with timescales of 2-6 years, 
caused by dumpiness of the accreting gas close to the black hole, is a likely explanation for the variations 
we see between the Chandra and XMM spectra. If this is indeed the case, and as already pointed out by 



Risaliti et al.l (l2002h . this scenario could be more common than previously thought, and applicable also to 
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intermediate Seyfert systems, such as Mrk 6. Observations of more such sources may lead to better con- 
straints on the geometry and composition of the gas distribution surrounding the AGN. Further observations 
at soft and hard energies, and detailed X-ray spectroscopy of this source would be of extremely useful for 
this purpose. 

The variable obscuration also poses some questions about the mechanism fueling the AGN. Our anal- 
ysis of the AGN core shows that the dominant contribution to the soft excess in this system comes from 
the AGN non-thermal emission, although there is some contribution from the jet, which could explain the 
flat power law spectra we are observing. Our results on the energetics of the source show that Bondi ac- 
cretion might not be enough to power it, so that it is possible that cold gas accretion is the dominant mech- 



anism in this system, u nhke what is seen in many more powerful radio sources (IHardcastle et all 12007 



Balmaverde et al 



2008h. Cold gas has also been suggested as the primary source of fuel for Cen A and 



NGC 3801 (|Croston et al.ll2007h . It is also possible that the active nucleus of Mrk 6 is accreting both the 
hot phase o f the ISM and co ld gas from a past merger. This is consistent with recent results on early-type 
galaxies by IPellegrinil (|2010l) . 



While the precessing jet scenario is challenging, due to the high torque required to rotate the axis 
of the black hole nearly 180° within a time of ~ 10^ years, required to create both perpendicular sets of 
bubbles, it cannot be ruled out. The dust lane hinted in the optical images could be a sign of a merger, 
which could have provided the necessary torque, and the evidence for obscuration variability near the black 
hole could be further evidence. It would be interesting to follow this hypothesis and find any correlations 
with other systems where a similar variation in absorption has been studied. However, this evidence can 
also be interpreted in terms of a cluttered and uneven environment, which can be inferred on bigger scales 
from the [OIII] emission distribution, causing the radio bubbles to expand asymmetrically and in diff'erent 
directions. Again, more evidence, both from future observations of Mrk 6 and from other sources, will be 
useful to achieve an accurate interpretation of the underlying physics behind the variable absorption and the 
asymmetries we have observed. 
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